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Abstract 

In this paper, we show that the coefficients of the E-characteristic polynomial 
of a tensor are orthonormal invariants of that tensor. When the dimension is 2, 
some simplified formulas of the E-characteristic polynomial are presented. A re- 
sultant formula for the constant term of the E-characteristic polynomial is given. 
We then study the set of tensors with infinitely many eigenpairs and the set of 
irregular tensors, and prove both the sets have codimension 2 as subvarieties in 
the projective space of tensors. This makes our perturbation method workable. By 
using the perturbation method and exploring the difference between E-eigenvalues 
and eigenpair equivalence classes, we present a simple formula for the coefficient 
of the leading term of the E-characteristic polynomial, when the dimension is 2. 
Keywords: E-Eigenvalues, tensors, E-characteristic polynomials, eigenpair equiv- 
alence class, irregularity 



1. Introduction 

Eigenvalues of higher order tensors were introduced in 2005 p___, [13] and 
have attracted much attention in the literature and found applications 
in science and engineering. The E-eigenvalues of a tensor were intro- 
duced in [T5l [17] . E-eigenvalues are invariant under orthonormal coordi- 
nate changes [17]. They have applications in determining positive defi- 
niteness of a multivariate system [Q13, IS]? best rank-one approximation 
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[T51 [TBI [T§1 [21] , magnetic resonance imaging [3, [2D] , spectral hypergraph 
theory [9] and symmetric embedding [21], and possess links with geome- 
try [H El 0, SI [16] . The numbers of E-eigenvalues and E-eigenvectors have 
been studied in [71, [Hj. For a real tensor, an E-eigenvalue with a real 
E-eigenvector is called a Z-eigenvalue. Z-eigenvalues are real and play the 
main role in the above applications. Algorithms for finding Z-eigenvalues 
were studied in [121 12] • 

An n-dimensional mth order hypermatrix A with entries in C is a map 
from {1, • • • , n} m — > C. We use notation A = (cLi 1 i 3 ...i m ) to denote such a 
hypermatrix. 

In physics [TT], [22] , a tensor is a physical quantity. In geometry, a tensor 
can be regarded as a multi-linear function in a Hilbert space, see discussion 
in Section 3. In such applications, in an orthonormal coordinate system, 
a tensor is expressed by a hypermatrix. It is important to determine 
parameters and properties of such a tensor, which are invariant under 
orthonormal coordinate changes. In Theorem 1 of Section of this paper, 
we will show that the coefficients of the E-characteristic polynomial of 
a tensor are orthonormal invariants of that tensor. For that discussion, 
we need to distinguish a tensor A and its hypermatrix A in a coordinate 
system. 

For some other applications, such as determining positive defmiteness 
of a multivariate system, best rank-one approximation, spectral hyper- 
graph theory and symmetric embedding, there are no coordinate systems 
involved. People just regard a hypermatrix A as a tensor. Hence, in the 
other parts of this paper, a tensor simply means a hypermatrix. 

For an mth order tensor A and a vector i G C", we define Ax m ~ l to 
be the vector in C n with its ith component as 

n 

{Ax )i = ^ ] a ii2-i m X i2 ' ' ' X i m - 

If A G C and x 6 C n satisfy 

(Ax m - l = \x, 

\ x T x = 1, U 
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then A is an E-eigenvalue of A and x is an E-eigenvector of A. If A and x 
are real, then A is also real [15]. In this case, A is called a Z-eigenvalue of 
A and x is called a Z-eigenvector of A. 

E-eigenvalues and Z-eigenvalues are invariant under orthonormal co- 
ordinate changes [T5]. An even-order real symmetric tensor always has 



Z-eigenvalues [15]. It is positive definite (semi-definite) if and only if 
all of its Z-eigenvalues are positive (nonnegative) [15] . The Z-eigenvalue 
with the largest absolute value and its corresponding Z-eigenvector form 
the best rank-one approximation of a real symmetric tensor [T5l [T9] . The 
modifier "E-" stands for Euclidean, as (CD) implies that the Euclidean norm 
of x is 1 if x is real. The modifier "Z-" names after Prof. Shuzi Zhou, 
who suggested (CD) to the author of [15]. Prof. Zhou departed by cancer 



in 2009. The name "Z-eigenvalue" has already been used widely in the 
literature 121 El SI El CZl IHl 021 031 HSl HBl QZZl HEl HSl EDI 1231 - Hence, 
we keep to use this name. 

E-characteristic polynomials were introduced in [TS, [17] , and discussed 
in [T4l [6] . Until now, there are no other papers containing discussion on E- 
characteristic polynomials. We now survey the definitions and properties 
of E-characteristic polynomials in P, dU, \T5[ [T7] 

Definition 1 Let A = (di^-i^) be an n- dimensional mth order tensor, if 
there exists x E C n \ {0} such that 



Ax m ~ l = 0, 
x T x = 0, 



(2) 



then A is called irregular. Otherwise, A is called regular flR \T_ 

When m is even, the E-characteristic polynomial iPa(X) is defined 175}/ 

as 

h(m,n) 

( m — 2 \ . 
Ax m - 1 -X(x T x) 2 x) = aj(m,n)\ J , (3) 

j=o 

where the second equality is the expansion of the resultant of in terms 
of X, h(m, n) is the highest power with generically non-zero coefficient and 
aj(m,n) 7 s are polynomials in entries of A. 
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When m is odd, the E- characteristic polynomial ^(A) is defined [F/\j 



as 

,m-l \™m-2„,\ h (m,n) 



^ A (A) = Res M ( M X \ = ]T a J (m,n)A 2 - ? , (4) 

\ x x x / 

where 2h(m, n) is the highest power with generically non-zero coefficient. 

Lemma 1 For a regular tensor, E- characteristic polynomial can be de- 
fined as 



Res 



f Ax m - 1 - \x^- 2 x 



for all m. 

Proof. For a regular tensor and even m, the systems 

f Ax™- 1 - \(x ) m ~ 2 x = 
\ x T x — Xq = 

and 

Ax 111 ' 1 - A(x T x)^x = 

are equivalent. □ 

According to Theorem 4 of [17] . any E-eigenvalue A of A must be a 
root of ^(A); if A is regular, then a complex number A is an E-eigenvalue 
of A if and only if it is a root of ^a(A). It was proved in [H] that when 
m is even, 

( 1 \ n 1 n—l 

Km, n) = [m ~ ] ~ = ]T(m - 1)\ (5) 



m 

i=0 



The definition (0Q) is not a strict extension of the classical definition 
for eigenvalues of a square matrix. The main point is that it excludes 
complex eigenvalues whose eigenvectors satisfy x T x = but x / [T7] , 
Recently, based upon this observation, Cartwright and Sturmfels [6] in- 
troduced equivalence classes of eigenpairs. If A G C and x G C n \ {0} 
satisfy 

Ax m ~ l = Ax, (6) 



then (A, x) is called an eigenpair of A. Two eigenpairs (A, x) and (A', x') 
are considered to be equivalent if there is a complex number t ^ such 
that A' = t m ~ 2 \ and x' = tx. Cartwright and Sturmfels proved that (Ej) is 
true for all m when counting the number of generalized eigenpairs without 
any normalization restriction [B] . 

An important property of a second-order tensor is that not only its 
eigenvalues but also the coefficients of its characteristic polynomial are 
invariants of that tensor. It was proved in [T5], [T7] that E-eigenvalues 
of a higher order tensor are invariants of that tensor. This implies that 
the coefficients of the E-characteristic polynomial, divided by the first 
nonzero coefficient, are invariants of that tensor, see Theorem 5 of |17j . 
There was no existing result about the invariance of the coefficients of the 
E-characteristic polynomial. 

In [6] , it was given that 

a (3,2) = (Res, (Ax 2 ))\ 

and mentioned in general that ao(m, n) is a power of Res, (^x"^ 1 ), with- 
out specifying the value of the power or giving a proof. When m = 3 and 
n = 2, they also present the coefficient of the leading term as the negative 
of the sum of two squares of linear combinations of the entries of A. 

In this paper, we explore more in this direction. The structure of the 
paper is as following. 

Since resultants are main tools to study the E-characteristic polyno- 
mials, we recall some facts about resultants in Section 2. In Section 
3, we prove that all the coefficients of the E-characteristic polynomial 
are orthonormal invariants of a tensor. We then show that when m is 
even, ao(m,n) is Kes x (Ax m ~ 1 ) , and when m is odd, it is the square of 
Res^Ax™- 1 ). 

In Section 4, we give two simplified determinantal formulas for the E- 
characteristic polynomial when n = 2. In the even case, this is a (2m — 
2) x (2m — 2) determinant. In the odd case, this is a (3m — 4) x (3m — 4) 
determinant. The parameter A only appears in the first m rows of the 
determinants. 



5 



We then address the leading coefficients of E-characteristic polynomi- 
als. Two sets of tensors are easy to handle: i) the set of tensors with 
infinitely many eigenpairs, ii) the set of regular tensors. Actually, the 
E-characteristic polynomials of tensors in the first set are identically zero; 
while a complex number is an E-eigenvalue of a regular tensor if and only 
if it is a root of the E-characteristic polynomial of that tensor. For ten- 
sors with finitely many eigenpairs, we apply the perturbation method, 
to approximate them by regular tensors with finitely many eigenpairs. 
Therefore we have to make sure that the perturbation is possible. Sec- 
tions 5 and 6 serve this purpose. In Section 5, we prove that the set of 
tensors with infinitely many eigenpairs has codimension at least 2 in the 
projective space of tensors. This shows it is possible to approximate a 
tensor in a hypersurface in the projective space of tensors by tensors with 
finitely many eigenpairs. We show in Section 6 that the set of irregular 
tensors has codimension at least 2 in the projective space of tensors too. 
Therefore, the union of the set of tensors with infinitely many eigenpairs 
and the set of irregular tensors has codimension at least 2 in the pro- 
jective space of tensors, which make the approximation to a tensor on a 
hypersurface by by regular tensors with finitely many eigenpairs possible. 

Finally, in Section 7, we introduce the deficit system to explore the 
difference of E-eigenvalues and eigenpair equivalence classes of Cartwright 
and Sturmfels [5]. Based on this analysis, we show that when n = 2, the 
leading coefficient of the E-characteristic polynomial is the ( II y^)th power 
of the sum of two squares when m is even, and the negative of the (m— 2)th 
power of such a sum when m is odd. 

2. Resultants 

First we recall some results for resultants since our discussion heavily 
depends on the properties of resultants, according to [S] and [ ITU] . 

For fixed positive integers di, • • • ,d n , and each pair (i, a), where i = 
1, • • • n, a = (ai, • • • , a n ) 6 Z> with \a\ = = ^ we introduce a 

variable Ui^ a . We have the following results about resultants. 
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Proposition 1 There is a unique polynomial Res x = Res^,...^ £ %[ui, a ] 
with the following properties: 

a. If Fi, • • • , F n G C[rci, • • • ,x n ] are homogeneous of degrees di, • • • , d n , 
then the equations F\ = • • • = F n = have a solution in CP" -1 if 
and only if Res{F\, • • • , F n ) = ; where the notation Res{F\, • • • , F n ) 
means Res(ci^ a ), if 



c. Res is irreducible in <C[ui )0e ]. 

Proposition 2 Res(Fi : ■ ■ ■ ,F n ) is a homogeneous polynomial in the co- 
efficients of Fi, with degree d\ • • • di-\di + \ ■ ■ ■ d n , for i = 1, • • • , n. 

Proposition 3 If Fj 7 s are homogeneous of degree d, and G{ = a ij^j) 
where (a^) is an invertible matrix, then 



Proposition 4 If Fj = F'-F" is a product of homogeneous polynomials, 



Res(F lr -- ,F jV -- : F n ) = Res(F 1: --- , F' p ■ ■ ■ , F n )Res(F h ■ ■ • , F'! ■ ■ ■ ,F n ). 



For an n-dimensional vector space V (over R. or C), let L : V — > V be a 
linear transformation, it induces a natural map L* on functions on V. Fix 
a basis {ei, • • • , e n } of V, Xi(^2 Cjef) = q defines a linear function X{ on V, 
for i = 1, • • • n. If the transformation matrix of L under this basis is (a^), 
then L*Xi = Y2 a ij x j- ^° f° r a polynomial F in variables it is 

a function on V, and L*F is a polynomial in variables x\, ■ • • x n obtained 
from F by substituting x^ with Y^ a ij x j- 

Corollary 1 For an invertible linear transformation L : V — > V, 




b. Res(x 



di 
l ' 



Res(Gi : ■ ■ ■ 



G n ) = det{a l3 ) dn ~ Res{F h --- ,F B ). 



Res(L*F ] 



L*F n ) = det(aij) 



Res(F h --- 



F n ) 
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Proof. Let 



G. 

Uj/yX . 



i, a 

a\=di 



Then Res(Fi : • • • , F n ) is an irreducible polynomial in variables Ui }0l , and 
by definition, if we denote 



L*Fi = ^ V i,a 



OC - 



then v^ a is a linear combination of u- variables Uj^, and Res(L*Fi, ■ ■ ■ , L*F n ) 
is obtained by substitute Mj jQ , by i;^ a , so it is a polynomial in variables 
and of the same degree as Res(Fi, ■ ■ ■ ,F n ). 

Since the matrix (a^) is invertible, the system Fi = ■ ■ • = F n = has 
a nontrivial solution if and only if the system L*F\ = • • • = L*F n = has a 
nontrivial solution, therefore Res(Fi, ■ • • , F n ) is a divisor of Res(L*Fi, • • • , L*F n 
We already know they have the same degree, so 

Res(L*F u • • • , L*F n ) = c Res(F u • • • , F„), 

where c is a constant. 

We can determine this constant by checking the case Fi = xf, i = 
1, • • • , n. By definition 

r* _ 

3 

Now by Proposition |U and the fact that Res\ y .. ,i is the standard determi- 
nant, 

Res{L*F u • • • , L*F n ) = Res((J2a ljXj ) d \ • • • , (Es^j) 4 ) 

xites((E>ijZ J -) dl - 1 J • • • , (E a n^j) dn ) 

= Res(Y,aijX j: (J2 a 2jX j ) d2 • • • , (E a ^) d ") rfl 

= Resi^aijXj, E a 2i^j, ■ ■ • , E a nj Xj) dv " dn 
= det(a i:j ) dr " dn . 
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Hence, 

c = det(a ij ) dv " dn . 

□ 



3. Orthonormal Invariance of the Coefficients 

The E-characteristic polynomial i/ja(^) is as defined in the introduction. 
Let a^( m n )(m, n) be the first nonzero coefficient of ^a(A). Then, for j = 
0, • • • , h(m, n) — 1, according to the relations between roots and coefficients 
of a one- variable polynomial, 

S h(m,n)-j - (- 1 ) /l(m ' n) ~ J aj(^^)/^(m,n)(^^) 

is the sum of all possible products of h(m, n) — j E-eigenvalues, hence an 
elementary symmetric function of the E-eigenvalues. 

In the following discussion, we distinguish a tensor A and its hyperma- 
trix A in a coordinate system. 

Let V be an n-dimensional vector space over R with a Euclidean inner 
product <, >. Let A be an mth order covariant tensor on V, that is, 

A : V m — ► R 

is a multi-linear function. In particular, a first order covariant tensor is 
called a covariant vector, denoted by x. In practice, what we consider is 
complexified tensor, i.e., to view them as multiple- linear functions from 

r%c^ c. 

We choose an orthonormal frame E = {ei, e n } in V : denote 

A(ei 1 , ei m ) .= a,i x i T ..i m , x(ej) = X{. 

Then (a^.-.^J (resp. Xi) is the hypermatrix of A (resp. x) with respect 
to the frame E. If all a^^...^, for ii, • • ■ , i m = 1, • • ■ , n, are real, then we 
say that A is real. If all Xi, for % = 1, • • • , n, are real, then we say that x 
is real. 
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If we choose another orthonormal frame E = {ei, e n }, that is 

where the matrix (C|) G 0(n) is the transformation matrix, then the 
hypermatrices of A and x with respect to ei, e n are given by 

^hii-im = 'y ■■■Cl m a j\j-2-j m i 0) 

Remark. For a tensor, the hypermatrices w.r.t. different orthonormal 
frames are orthogonally similar. 

Remark. Since we only consider tensors under orthonormal transforma- 
tions, we do not distinguish covariant tensors and contravariant tensors. 

For a complex mth order covariant tensor A and a covariant vector x, 
we can define a covariant vector Ax. m ~ l by 

n 

Ax m -\ei) = A ^ e ^ ' • • ' e ^>(e* 2 ) • • • x(e im ), (9) 

for i = 1, • • • ,n, and any orthonormal frame £" = {ei, e n } in V. For 
any covariant vector x, define 

n 

x r x = ^(x( ei )) 2 , (10) 

i=l 

for any orthonormal frame E = {ei, ...,e n } in V. By (J7j) and (jSJ), we see 
that the definitions (EJ) and ( TTOT ) are well-defined, independent from the 
frame E. 

For a complex mth order covariant tensor A , if A G C and non-zero 
covariant vector x satisfy 

f Ax m 1 = Ax, 
I x^x=l, 

then A is an E-eigenvalue of A and x is an E-eigenvector of A. By the 
above discussion, we see that A is invariant, i.e., it is independent of the 
choice of the frames. 
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We say that a fact or quantity is orthonormal invariant, if it is 

invariant under the changes of orthonormal frames. 

Fix an orthonormal frame E = {ei, ...,e n } in V, let the hypermatrix 
A = (a^ - im) ( res P- x i) be the hypermatrix of A (resp. x), let 

n 

il"- ,Vi=l 

we use the notation Res x (*4x m_1 ) to denote Res, (Fi, F n ), and ^/u(A) 
to denote ^a(A). 

Theorem 1 Res x (ix™ 4 ) is an orthonormal invariant of A. Further- 
more, all the coefficients of ipj[(X) are orthonormal invariants of A. 

Proof. We first prove that Res x (ix™ 4 ) is an orthonormal invariant 
of A. 

If we choose another orthonormal frame ei, e n , 
and denote 

n 

Fi(x) = ^ ] CLa 2 ...i m Xi 2 • • ■ Xi m , 

12,— ,irra = l 

by a direct calculation we get 
So if we let 

G l (x) = Y J ClF ] (x), 

then by Proposition El 

Res, (Gi, G n ) = detiCl)^- 1 ^ 1 Res x (F h F n ) , 
and noticing 

Fi(x) = Gi(^2 Cjxj, • • • , ^2 CjXj), 
j j 
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by Corollary of Proposition H[ we have 

Re Si (A, F n ) = <M(Cf ^"^Res, G n ) . 

Since det(Cj) = ±1, (m — is even, we have 

Resj (P u Fnj = Res x (F b F„) . 

By the same method we can prove that when m is even, Res x (Ax" 1- " 1 
—X (x xj 2 x) is orthonormal invariant, when m is odd, Res( X)Xo ) I x r x _J 
is orthonormal invariant. The proof is completed. □ 

Because of this theorem, in further discussion, we do not distinguish a 
tensor and its corresponding hypermatrix strictly. 

Lemma 2 For a hypermatrix A, Res x (Ax™ -1 ) is an irreducible polyno- 
mial in the entries of A (view them as variables). 

Proof. For any index a = (ai, • • • , a n ) with |a| = m — 1, it defines a 
subset S a of {1, • • • , n} TO_1 , 

S a = {(^2, • • • , im) |among i^ ■ ■ ■ , i mi ctj of them are equal to j, for j = 1, • • ■ , n}. 
Obviously, 

S a Pi S/3 = 0, for any a^, 
U a S a = {l,--- ,n} m -\ 
For the system Ax m ~ l = 0, if we denote the equations by 

i = 1, • • • , n, then 

Assume Res x (^x TO_1 ) = P(q a ) = f(a il ...i m )g(a lv .. im ). Let us fix a 
function 

(/) : = (ai, • • • , a n ), |a| = m — 1} — >• {1, • • • , n} TO_ , 
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such that 4>(a) G S a . We can choose entries of A such that for every Q a 
only the one entry nonzero is aw a ). Then Kes x (Ax 171 " 1 ) is an irreducible 
polynomial in the variables corresponding to these nonzero entries. Hence, 

P{^i<t>{a)) = /(O»0(a))p(fli0(a))- 

Since P is irreducible, one of /(a^a)), ^(^(a)), say <7(a^( a )) has to be 
of degree and the degree of the other equals to the degree of P. Then 
g(ai v ..i m ) has degree because P(a,i r ..i m ) has the same degree as P(ci a ). 
Therefore Res x (Ax 771 - 1 ) is irreducible. □ 

We now give a general formula for ao(m, n). 
Theorem 2 When m is even, we have 

a {m,n) = Res x (Ax m - 1 ) . (11) 
When m is odd, we have 

a (m,n) = (Res x (Ax™- 1 )) 2 . (12) 
Proof. When m is even, by (E3), we have 

a (m, n) = ^(0) = Res x (Ax m ~ l ) . 

This proves (TTTT) . 

When m is odd, by (ED, we have 



a (m, n) = i/j a (0) = Res^) T ^ 



To prove (JT21). it suffices now to prove 



/ A r m-1 \ 

Hes^o) t 2 = (Res,^- 1 )) 2 . (13) 



We first show that (Res^ (Ax 171 - 1 )) has the correct degree. Denote 
(Ax 777 - 1 ), as Fi(x) for z = 1 • • • n. Since x T x — x\ has degree 2, by 

Proposition El for any i, ReS( X;Xo ) fjr^^J is a homogeneous polynomial 
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in the coefficients of Fi(x), with degree 2(m — l) n 1 . Obviously, for any i, 
(Res^ (Ax™- 1 )) is a homogeneous polynomial in the coefficients of Fi(x) : 

with degree 2(m — This shows that (Res^ (Ax 171 ' 1 )) has the correct 
degree. 

Next we prove that the system 

has a nonzero solution if and only if Res^ (Ax 771 ^ 1 ) = 0. Let (x, Xq) be a 
nonzero solution of flHJ). Since A is regular, xq ^ 0. By the last equation 
of ( 1141) . x 7^ 0. Then £ is a nonzero solution of 

Ax 771 ' 1 = 0. (15) 

Thus, Res x (Ax 771 ^ 1 ) = 0. On the other hand, suppose that Res x (Ax™' 1 ) 
= 0. Then (1151 ) has a nonzero solution x. Let Then (x, xq) is a 

nonzero solution of (HI). Therefore Res M ( = 0andRes x (Ax 771 - 1 ) 

= define the same varieties. By the irreducibility of Res x (Ax 771 ^ 1 ) , we 
know 

ReS (-^o) ( JL_ , J = ( ReS - (^ m_1 )) fc - 

Jb Jb Jb^J 

Then by degree, we know k = 2. These prove (1T31). The proof is completed. 
□ 

This shows the importance of Res x (Ax 771 ' 1 ). When m = 2, it is the 
determinant of square matrix A in the classical sense. Hence, it is a 
genuine extension of the determinant of a square matrix, and deserves to 
be studied further. 

Before ending this section, we give a proposition on the degree of the 
coefficients of ^(A) as polynomials in the entries of A. 
When m is even, we may see that 

m-2 

Ax m - 1 -X(x T x) 2 x = (16) 
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is a system of homogeneous polynomials in x. Every equation of ( TTBT ) has 
the same degree m — 1. Thus, by Proposition El Vu(A) is a homogeneous 
polynomial in the entries of A and A, with degree n(m — l) n . Hence, in 
(EJ), aj(m, n) is a homogeneous polynomial in the entries of A, with degree 
n(m — l) n_1 — j. When m is odd, the first m equations of 

( Arr-m-l _ \ m-2 

\ T 2 17 

have the same degree m — 1, while the coefficients of last equation of 
( JT7I ) are either 1 or —1. Thus, by Proposition [21 ^a(A) is a homogeneous 
polynomial in the entries of A and A, with degree 2n(m — l) n . Hence, in 
(Sj), a,j(m, n) is a homogeneous polynomial in the entries of A, with degree 
2n(m — l) n ~ l — 2j. We now have the following proposition. 

Proposition 5 In (Ejj, a 7 (m, n) is a homogeneous polynomial in the en- 
tries of A, with degree n(m — l) n ~ l — j. In particular, dh(m,n){' r ^i n ) ^ s a 
homogeneous polynomial in the entries of A, with degree n(m — 1) — 

(m-l)"-l 

771-2 

In (2J) ; aj(m,ri) is a homogeneous polynomial in the entries of A, with 
degree 2n(m — l) n — 2j . In particular, ah(m,ri){ m i n ) ^ s a homogeneous 
polynomial in the entries of A, with degree 2n(m — l) n ~ l — 2 ^ m m ^ 2 — 



4. E-Characteristic Polynomials when n = 2 

We now derive some simplified forms for ^(A) in the case that n = 2. 

Let bj = J2{ a U2-i m '■ exactly j - 1 of i 2 , • • • ,i m are 2} for j = 1, • • ■ , m, 
and Cj = YsW^T-hn '■ exactly j - 1 of i 2 , • • • ,i m are 2} for j = 1, • • • , m. 
For j = 1, • • ■ , m — 1, let = &j — Cj+i, and let d m = 6 m . 

For an N x TV matrix M = (m^), just like {ran,m22,-"" , ^iwv} is 
its diagonal, we call {mi^fc+i, ^2,2fc+2, ■ • • , ^7V-2fc,7v} the /cth even upper 
sub-diagonal of M for positive A: such that 2k + 1 < N. 

We now discuss the case that m is even. 
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Theorem 3 Suppose that m = 2k + 2 and n = 2, where k > 0. Lei 

&2j+l = 02j+l ~ 



yj A and c 2j +2 = c 2j+2 - Q j A, 



/or j = 0, • • • , k. Then for a regular A, ^(A) is the determinant of the 
following (2m — 2) x (2m — 2) matrix: 



Mi = 



bi 


b 2 


6 3 •• 


• b m -\ 


b m 





•• 


■ 





h 


62 •• 


■ b m -2 


bm—1 


bm 


•• 


• 








•• 


■ h 


b 2 


h 


64 •• 


bm 








•• 


ci 


C2 


C3 


c 4 • • 


Cm 


-ci 


o!i 


d 2 ■ ■ 


■ d m -2 


djji — 1 




•• 


■ 





-ci 


d\ ■ ■ 


dfn—3 


d m -2 


d<m—l 


bm 


■ 








•• 


d\ 


d 2 


d-s 


C?4 • • 


bm 



(18) 



m-2 



7 



even upper 



Here, the first m — 1 entries of the diagonal and the first 
sub-diagonals of Mi are 61, 63, • • ■ , 6 m -i, tuMe £/ie mf/i entries of the diag- 
onal and the first even upper sub-diagonals of Mi are c 2 , C4, • • ■ , c r 
They are linear factors of X. 

Proof. Consider 



5 ' 



m^2 
2 



F(x,A) = 



(^x m - 1 ) 1 - A (x T x) 2 xi 
x 2 (Ac™- 1 ^ - a?i (Ac™- 1 ), 



= 0. 



(19) 



Since A is regular, any nonzero solution of ( TTBl ) is a nonzero solution of 
(EE5D- This implies that ^a(A) is a factor of Res x F(x, X). On the other 
hand, the only possible additional nonzero solution of (JTSj) satisfies Xi = 
and (^x m - 1 ) 1 = 0. Since (Ac™" 1 ^ = E£i & i<~*4 _1 > £1 = 0, x 2 ^ 
and (Ar 1 ™ -1 ^ = imply that b m = clyi-2 = 0. Therefore, we conclude 
that 

Res x F(x,A) = 6 m ^(A). (20) 
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Now, we have 



F(x,\) 



Ek m—2j—l 2j n 

7=0 2j + lX 1 ■ X 2 + 02j+2X 1 



„m-2j-2 2j+l 



X, 



By the Sylvester formula P [TO], Res x F(x, A) is the determinant of the 
following (2m — 1) x (2m — 1) matrix: 



bi 


&2 


63 •• 


• b m -\ 







•■ 











&1 


62 •• 


■ b m - 2 


bm-1 


bm 


•• 














•• 


■ h 


b 2 


h 


64 •■ 


bm 











•• 





h 


b 2 


63 •■ 


" b m -l 




-ci 


rfl 


d 2 ■ ■ 


■ d m -2 


d m -l 


b m 


•■ 











-ci 


d\ ■ ■ 


dm— 3 


d m -2 


dm—1 


bm 














•• 


d\ 


d 2 


d 3 


G?4 • ■ 


bm 











•• 


■ -ci 


di 


d 2 


d 3 ■ ■ 


dm—1 


bm 



As elementary row and column operations preserve the value of the de- 
terminant, we may subtract the last row from the mth row in the above 
matrix. Then all the elements of the last column of the matrix are zero 
except the bottom element is b m - Then we may delete the last column 
and the last row, and extract b m from the determinant. Then we see 
that Res x F(x, A) = 6 TO det(Mi). Comparing with (l2~Uj). we have ^a(A) = 
det(Mi). The proof is completed. □ 

The merit of M\ is that only the first m entries of the diagonal and the 
first ' nk Y L even upper sub-diagonals of M\ contain linear factors of A, i.e., 
those bi and q. By some elementary column operations, we may eliminate 
A in the terms other than b\ and c-2- 

When m = 4 and n = 2, we have /i(4, 2) = 4, 



^(4,2)^4(4, 2) = (61 - c 2 - 63 + c 4 ) 2 + (ci + 6 2 - c 3 - 6 4 ) 2 
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+ ^2222 — (21221 — «1212 ~ &1122 _ «2211 ~ ^2121 _ 02112] 



+ {ai211 + &1121 + ^1112 + ^2111 — &1222 — «2221 _ ^2212 _ ^2122) • 

We now discuss the case that m is odd. Let e\ = b\C\,e 2 = b\c 2 + 
hPu e 3 = hcs + &2C2 + &3C1, • • • ,e m = bic m + 6 2 c m -i H h b m c u e m+ i = 

b 2 C m + • • • + b m C2, • • • , e2 m -l = & TO C m . 

Theorem 4 Suppose that m = k + 2 and n = 2, where k > 1 is odd. Let 



e 2j = e 2j 



X 



.J 



for j = 1, ■ • ■ , m - 1. Let ft = e l + 6i^+i /or i = 2, 4, • • • , m - 1; 
ft = ei + hd l+1 for i_= 3,5, • • • ,m - 2; ft = e { - c m d;_ m +i for % = 
m, m + 2, • • • , 2m — 3; fi = e~i — c m di- m+ i for i = m + 1, m + 3, • • • , 2m — 2. 
TTien «s t/ie determinant of the following (3m — 4) x (3m — 4) matrix: 

\ 



M 2 = 



h 


h ■■ 


fm—l 




• e2 m -2 


G 2m -i ■ • 





ei 


e 2 • ■ 


■ e m _2 


6m— 1 


• e2m-3 


&2m-2 ■ ■ 








•• 


e 2 


e 3 ■ ■ 


fm—1 


fm 


' /2m-2 


-ci 


d\ 


dm—l 


bm 














-Ci ■ ■ 


■ d m -2 


dm—l 














•• 





•• 


d 2 


d 3 








•• 





•• 


d\ 


d 2 


b m 



) 



21) 

Here, the first m — 2 entries of the first row and the last m — 2 entries of 
the rath row are somewhat different from the other entries of the first m 
rows, with or e~i being replaced by fi or fi . 



Proof. Consider 



G{x,\) 



x 



x 2 (Ax™- 1 



% - A 2 
) 1 - xi (Ax m - V 



x T x) 



2m-4 



X\X 2 



= 0. 



'22^ 
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Similar to the proof of Theorem EH we may conclude that 



Res x G(x, A) = b m ci^ A {\)- 



(23) 



Now, we have 



GO, A) = 



E 



m-l 2to-2j'-2 2j , - 2m-2j-l 2j— 1 ' 

j=Q e 2j+l x \ X 2 + Z^j=l e 2j^l ^2 



By the Sylvester formula P HD], Res x G(x, A) is the determinant of the 
following (3m — 2) x (3m — 2) matrix: 

/ 



ei 


G2 


e 3 • ■ 




6m 


• e2m-2 


62m-l - " 











ei 


e 2 • ■ 


• e m _2 


6m— 1 


• e2m-3 


e~2m-2 • • 














•■ 


e 2 


e 3 • ■ 


6m— 1 




• e~2 m _2 


62m- 1 


-ci 


d\ 


d 2 ■ ■ 




■■ 

















-ci 


d\ 


djn—\ 






















-ci • ■ 


' <^to-2 


dyn—1 




















•■ 





•■ 




d 3 














•• 





•■ 


di 


d 2 


b m 











•■ 





•■ 


• -ci 


d\ 


dm—l 


bm 



y ••• —c\ d\ • • • a TO -i b m J 

Note that e\ = &1C1 and e2 TO -i = b m c m . Thus, as elementary row and 
column operations preserve the value of the determinant, we may eliminate 
ei in the first column with a multiple of the (m+l)th row, and eliminate 
e2 m -i in the last column with a multiple of the last row. Then we may 
delete the first column, the last column, the (m + l)th row and the last 
row, and extract b m c\ from the determinant. We see that ~Res x G{x, A) = 
6 TO cidet(M2). Comparing with (1231). we have ^a(A) = det(M2). The proof 
is completed. □ 

The merit of M2 is that only the first m entries of the diagonal and 
the first m — 2 even upper sub-diagonals of M2 contain linear factors of 
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A 2 , i.e., those &{ and fa. By some elementary column operations, we may 
eliminate A in the terms other than fa and e 2 . 
When m = 3 and n = 2, in [6], it was given that 

&/i(3,2)=3(3, 2) = — (—am + ai 22 + a 2i2 + a22i) 2 — (am + &121 + «2ii — fl222) 2 

If we use 6^ and q to write it, then we have 

a/*(3,2)=3(3, 2) = -(&!- c 2 - 6 3 ) 2 - (ci + 6 2 - c 3 ) 2 . 

We see that in the form using bi and q, a^^) (3, 2) and (4, 2) are very 
similar. Let P m and Q m be the sum of the first m terms of the following 
two series respectively: 

bi - c 2 - h + C4 + h - c 6 - b 7 + c 8 + b 9 - c 1Q 

and 

ci + b 2 - c 3 - 6 4 + c 5 + 6 6 - c 7 - b 8 + c 9 + 610 , 

where the signs of these terms are cyclically changed with the cycle four. 
Then, we may write 

a M3 , 2) (3,2) = -(P 3 2 + Q 2 ) 

and 

a M 4, 2) (4,2) = P 2 + Q 2 . 
Using Theorem EH, we find via computation that 

a M6 , 2) (6,2) = (P 2 + Q 2 ) 2 . 

It seems that there are formulas for CLh(m,2){ m i 2) via P^+Q^ for all m > 2. 
Later we will show that this is true. 

5. A Bound on the Codimension of the Variety of 
Tensors with Infinitely Many Eigenpairs 

For positive integers di, g^, • • • ? dm let P be the complex projective space 
corresponding to homogeneous coordinate ring 

C[u ia I \a\ = a H h a n = d i: i = 1, • • • n}. 



20 



Now for a system of homogeneous polynomial equations: F\ = F2 = 
• • • = F n = where Fi, F2, • • • , F n are homogeneous polynomials of posi- 
tive degrees d\, d^ ■ ■ ■ , d n in variables 



Then modulo obvious scalar multiplication of systems, there is one-to-one 
correspondence between such systems and points in P. 

For such a system, the solution set in CP" has dimension > n — n = 0. 
For systems corresponding to generic points in P, the solution sets have 
dimension 0, i.e. they have only finitely many solutions in CP". 

Let X C P be the set of points whose corresponding systems have 
infinitely many solutions in CP", i.e., 

X = {([q j( J) I F\ = • ■ ■ = F n = has infinitely many solutions in CP"}. 

Obviously codimX > 1 in P, where X is the Zariski closure of X. 

Let us analysis the subset X in more detail. We will use a kind of 
hidden variable argument similar to [8], and we need some careful treat- 
ment to apply results on resultants, because the systems should have n 
homogeneous equations in n variables or n inhomogeneous equations in 
n — 1 variables. 

Fix a general system F\ = F2 = • • • = F n = 0, 



where a = (ao, cki, • • • , a n ). 

For < s 7^ t < n, we introduce a set of variables yj, j = 0, • • • s, • • • n 
(more precise notation would be Then in Fi, % = 1, • • • , n, we replace 
Xj by tjj for j 7^ s, and replace x s by ^rVt- Thus write Fi as a homogeneous 
polynomial of degree di in yj, j ^ s with coefficients in C[^f]: 
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(5 = (/3o, • • • , f3 s , ■ ■ ■ , f3 n ), and we denote it by F/'*. The set of multi-indices 
a for x- variables with a s = is one-to-one corresponding to multi-indices 
(3 for ^-variables. Obviously 

Let us make this treatment clear by an example. For the homogeneous 
polynomial F = x\x\ + xqx\ + xiic 2 ; + a; 2 ^, if we take s = 0, t = 1, we 
obtain the homogeneous polynomial [(f^) 2 + f^]2/i + 2/12/2 + 2/i2/2 after the 
replacement. 

The system Fj 8 ' = • • • = F*'* = is a system with n homogeneous 
equations in n-variables, so the resultant Res s,t = Res(F^ ,t , • • • , F*'*) is a 
polynomial in variables f*'p, thus a polynomial in Z[itj iQ ,, i = 1, • • • n, ^]. 

Lemma 3 J/^/ie system F\ = F 2 = • • • = F n = /ms a solution (ag, fli, ■ " ? 
wzt/i <2i 7^ ; then for any s ^ t, — is a solution o/i?es(F 1 s ' i , • • • , F*'*) = 0. 

Proof. The system J2\ a \ =d . ftt(%)y a = 0, i = 1, • • ■ , n has a non-trivial 
solution (fj)# fl , so J is a solution of Res(F*'\ • • • , F*>*) = 0. □ 

Lemma 4 For any (cj )Q! ) G X ; a£ least one of the polynomial Res s,t (ci ja , ^) 
vanished identically. 

Proof. Otherwise, a solution (a\,--- ,a n ) of the syetem F\ = F2 = 
• ■ • = F n = must satisfy Res s,t (ci ja , j-) = for any £ with 7^ and 
s 7^ t, such solutions can have at most finitely many. □ 

Theorem 5 The Zariski closure X of X has codimension at least 2 in¥ 
ifn>2. 

Proof. For < s ^ t < n, let 

*''*= j[ci >a ] lite* - '* ( c ^^) = 
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obviously X C UX S ' , so if we can prove that every X s ' has codimension 
at least 2, then so is X. 

As a polynomial in variable ^ with coefficients in Z[i^ )Q! , i = 1, • • • n], 
the constant term of Res s,t (ui^ a , |*) is the resultant for the reduced system 
of our general system by letting x s = 0, which is a general system in n 
variables yj, j ^ s. We denote the constant term by Res{ui^ a ^a s = 0), by 
formula (l24l) . it is an irreducible polynomial in variables ui a (with a s = 0), 

Now let us consider the coefficient of linear term in Res s,t (ui ia ,^r). 
Notice by formula (T2~4l), if we take the degree ^ to be 0, then Res s,t (ui ja , ^) 
is a homogeneous polynomial. 

For the coefficient of linear term, again by formula ( 1241) . every monomial 
is a monomial in variables Ui j0l with a s = times a variable ui^ a where 
a s = 1, so it has lower degree in variables u^ a with a s = 0, and it has 
the same degree as the constant term, so it is coprime with constant term 
(constant term is irreducible). 

Inside P (dim(¥) > 2 if n > 2), the subvariety defined by the constant 
term and the coefficient of linear term cannot be the zero set of a single 
polynomial, so must have codimension > 2. The subvariety X s,t is con- 
tained in this variety, so X s,t has codimension at least 2. □ 

Up to a common scalar multiplication, the set of hypermatrices is one- 
to-one corresponding to points in CP n ™ -1 . 

Theorem 6 In the space CP n ~ l , the set of hypermatrices with infinitely 
many eigenpairs has codimension at least 2 (n > 2). 

Proof. Let A be an mth order hypermatrix. Then by a similar argu- 
ment as in Theorem El the subset of hypermatrices A for which Ax m ~ l = 
has infinitely many solutions is contained in a subvariety (we denote it by 
Xq) which is defined by two coprime polynomials, and Xq has codimension 
> 2, i.e., the closure of the subset for which is an eigenvalue in infinitely 
many eigenpairs only is a codimension > 2 subvariety. 
Now for the system 

Ax == Ax, 
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it has infinitely many solutions with A 7^ if and only if 



Ary.m-1 _ to-2 

JxJu — 

has infinitely many solutions with xq 7^ 0. 

For s = 1, • • • n, we introduce variable y = 0, • • • , y s , • • • , y n , and 
rewrite the system Ax m ~ l = x™~ 2 x: replace Xj by yj for j 7^ s, and 
replace x s by ^yo, thus we have a new system Ff = F 2 S = • ■ ■ = F* = 0, 
this system is obtained from Ax m ~ l = by a replacement Xj by yj for 
j 7^ s and x s 1— > ^2/0 an d a translation of coefficients: with an extra term 
—y™~ 2 yi for i ^ s and — fj2/o" 1 f° r * = s - 

The resultant Res s (a iv ..i n: of the system F* = F| = • • • = = 
is a polynomial in Zfa^...^, let 




Then by similar argument in the proof of theorem [3, the set of ^4 for which 
Ax 711 ' 1 has infinitely many solutions with xq 7^ is a subset of 

Ug =1 X s , and X s has codimension at least 2. 

Thus, the set of hypermatrices with infinitely many eigenpairs is a 
subset of the subvariety Xq U" =1 X s . Hence its closure has codimension 
> 2. □ 



6. A Bound on the Codimension of the Variety of 
Irregular Tensors 

Let 

X = {[A]\A is irregular} C CP"™ -1 , 
Theorem 7 The Zariski closure X of X has codimension at least 2. 
Proof. For a tensor A, let 

Fi = ^ ] aa 2 ...i m Xi 2 ■ ■ • X{ m . 
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Then A is irregular if and only if 

F 1 = F 2 = ■ ■ ■ = F n = x T x = 

has a solution in CP 71 , this is a system with n+1 homogeneous equations 
in n variables. 

Now assume A is irregular, then for any i = 1, • • ■ , n, the system 

F x = F 2 = ■ ■ ■ = F t = ■ ■ ■ = F n = x T x = 

has a nontrivial solution, therefore the resultant Aj of this system is 0. 
For a hypermatrix A with only nontrivial entries OLjj-jj = 1, the system 

F x = F 2 = ■ ■ ■ = F t = ■ ■ ■ = F n = x T x = 

is 

\ x T x = 

it has no solution in CP n_1 , so Aj ^ 0. 

Notice for the resultant Aj, the variables aa 2 ...i n are missing, so the 
g.c.d of Ai, • • ■ , A n is 1. i.e, any component of the subvariety of CF 12 "^ 1 
defined by 

Ai = • • • = A n = 

cannot be defined by a single polynomial. Hence any component of the 
subvariety Ai = • • • = A„ = has codimension > 2. 
Now we have the conclusion because 

Xc{A 1 = --- = A n = 0}. 

□ 

7. The Leading Coefficient 

We now study the properties of the leading coefficient ah(m,n){ m i n ) by 
exploring the difference of the definitions of E-eigenvalues and eigenpair 
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equivalence classes. We see that beside those eigenpair equivalence classes 
which correspond to E-eigenvalues, all the other eigenpair equivalence 
classes correspond to nonzero solutions of 

/ Ax"*" 1 = Ax, 

1 X T X = , ( 25 ) 

We call (J25D the deficit system of A. 

Now we investigate the leading coefficient for the characteristic poly- 
nomials for tensors when n = 2. 

Lemma 5 For a regular tensor, the leading coefficient of its characteristic 
polynomial is 0, if and only if the system 

f Ax m - 1 = Xx 
\ x T x = 

has a nontrivial solution. 

Proof. Consider the inhomogeneous system 

f Ax™' 1 = Xx, 
1 x T x = 1. 

This system has no solution at oo because of regularity. 

If we consider A as a constant and homogenize the system with respect 
to xq, we have the following system 

f Ax m - 1 - Xx^- 2 x = 0, 
1 x T x — Xq = 0. 

The resultant for this homogeneous system (the E-characteristic polyno- 
mial) has the exact information of multiplicity of eigenpairs, i.e., it is a 
constant times [8] 

n (A ^ } (A-A0™^, 

where the product is over distinguish eigenpairs {(A^,^)}, and m(Xi,Xi) 
is its multiplicity. 
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Therefore, if the leading term is 0, there is some eigenpair missing, 
which has to be a solution of 

f Ax m - 1 = Xx, 
\ x T x = 0. 

On the other hand, if the leading term is not 0, the deficit system 
cannot have a nontrivial solution, otherwise A will have more eigenpairs 
than expected. □ 

Lemma 6 For a regular tensor, the system 

f Ax 712 - 1 = Xx, 
\ x T x = 

has a nontrivial solution if and only if 

where P m and Q m are the sum of the first m terms of the following two 
series respectively: 

h - c 2 - h + c 4 + h - c 6 - b 7 + c 8 + b 9 - c 10 

and 

Ci + b 2 ~ C 3 - 64 + C 5 + 6 6 - C 7 - 6 8 + Cg + b W , 

where the signs of these terms are cyclically changed with the cycle four, 
and bj = ^2{a lt2 ... im : exactly j - 1 ofi 2 ,--- ,i m are 2} for j = 1, • • • , m, 
and Cj = J2{ a 2i 2 -i m '■ exactly j -1 of% 2 ,--- ,i m are 2} for j = 1, • • • , m. 

Proof. Since n = 2, if the deficit system ( 1251 ) has a nontrivial solution, 
then all the nontrivial solutions of the deficit system ( l2~ol ) are nonzero 
multiple of (1, >/— T) or (1, — y/— 1). Therefore, there are at most two 
eigenpair equivalence classes corresponding to nonzero solutions of the 
deficit system ([25]), when n = 2. Note that when n = 2, Ax" 1 ^ 1 = Xx can 
be written as 

( h rrm-i i-l \ 
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Substituting (1, to pBj), we have 



tzm^ = a, 

Eliminating A (A ^ by regularity) , we have 

TO TO 

^(.,:(V Z T)"+E C *(^ T )*= a 

i=l i=l 

By the definitions of P m and Q m , we have P m = — Q m V~~^- Similarly, 
substituting (1, to ( 1261 ) and eliminating A, we have P m = Q m \f— 1. 

These imply + Q m = 0. 

On the other hand, if P^ + Q 2 = 0, by a similar but reverse argument, 



we see that (1, or (1, — y/— 1) is a nontrivial solution of the deficit 

system (1251). This completes the proof. □ 

Theorem 8 If A has only finitely many equivalence classes of eigenpairs, 
then the leading coefficient a^ m 2){m^ 2) = if and only if 

P 2 + n 2 n 

* m ' ^to u- 



Proof. First, if A has only finitely many equivalence classes of eigen- 
pairs, then its characteristic polynomial cannot be identically 0. 

Let us consider the subvarieties of CP 2 ™" 1 : X = {au m ,2)ijn : 2) = 0}, 
and Y = {P^ + Q 2 m = 0}, they both have codimension 1. Let Z\ = 
{[A]\ A has infinitely many eigenpairs}, and Z<± = {[A]\ A is irregular}. 
Then both Z\ and Z^ have codimension > 2. So X— Z\ — Z^ and Y — Z\ — Z2 
are not empty. 

By Lemmas [5] and E| X — Z\ — Z^ C Y and Y — Z\ — Z2 C X, taking 
closure, we have X C Y and Y C X. This implies X = Y and the 
conclusion follows. □ 

Theorem 9 If A has only finitely many equivalence classes of eigenpairs, 
then the leading coefficient a<h{m,2)( m -> 2) is 



Qm) 2 



m— 2 
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for even m, and 

-(pl+Qir- 2 

for odd m. 

Proof. Denote that 

W = {[B]\I»+<$ n = 0}. 
W can be written as a union of two irreducible varieties 

W = w l {Jw 2 , 

where 

W 1 = {[B}\P m + V^lQ m = 0}, 
W 2 = {[B}\P m -V^lQ m = 0}. 
Since P m + \f—lQ m and P m — \/—lQ m are irreducible, by Theorem |S] 

a h {m,2){m, 2) = (P m + >/=lQ m )*(P m - v^lQ™)^, 

where C is a constant. Since au m ^) (m, 2) G Z[a« r ..i m ], /c must be equal 
to s. Then by checking the degree, we see that our results hold up to a 
scaling constant. By checking an example such as a\...\ = a2- -2 = 1 but all 
the other entries of A are zero, we have the results. □ 

Acknowledgment We are thankful to two referees, whose comments 
greatly improved our paper. 
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